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Abstract 

This Article is a brief review of coarsening phenomena occurring in systems where quenched features - such as random 
held, varying coupling constants or lattice vacancies - spoil homogeneity. We discuss the current understanding of the 
problem in ferromagnetic systems with a non-conserved scalar order parameter by focusing primarily on the form of 
the growth-law of the ordered domains and on the scaling properties. 
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1. Introduction 


Coarsening can exhibit rather different features in systems where homogeneity holds or not. In order to keep 
the discussion at the simplest level I will focus on the class of ferromagnetic systems with a scalar order-parameter 
(spin), such as those whose properties can be described by the Ising model or, in a continuum approach, by the usual 
Ginzburg-Landau free energy. Inhomogeneous systems with a vectorial order-parameter are much less studied and 
will not be discussed here. Besides, the attention will be restricted to a dynamics that does not conserve the order- 
parameter. This descriptions are in principle suited, e.g., for a magnetic solid. The aim of the article is not to provide a 
comprehensive description of the behavior of existing real systems, but rather to discuss at the simplest possible level 
the new features introduced by the presence of inhomogeneities. 

In the kind of systems we are considering coarsening is usually observed after a quench from an high-temperature 
equilibrium configuration to a temperature below the critical one. Relaxation toward the new equilibrium state then 
occurs by the formation and growth of domains inside which the system is basically equilibrated in one of the two 
symmetry-related low-temperature equilibrium phases (see Article by L.F. Cugliandolo in this Volume). For a magnet, 
these are the two possible ordered configurations with positive or negative magnetization ±m. The non-equilibrium 
behavior occurs on the interfaces which move in order to increase the typical domain’s size. A pictorial representation 
of the process is depicted in hgureQ] 

The distinguishing characteristic of coarsening in homogeneous systems is the scaling symmetry 01 ■ Accordingly, 
in the late stage of the process, configurations of the system at different times are statistically equivalent if lengths are 
measured in units of the characteristic size L(t) of the domains, which grows algebraically like 

L(t) = at l,z , (1) 

where a is a constant, and z — 2 in the case without order-parameter conservation considered here |T|]. The scaling 
symmetry is mirrored by the form taken by different quantities, as for instance correlation functions. Using the 
language of spin systems the two-time spin-spin correlation G{r, t) = (S i(t)S where S , = ± 1 are spin variables 
on sites i and j of a regular lattice at a distance r, obeys 


G(r,t w ,t ) = g 


r Lit) 
L(t w ) L(t n ) 


( 2 ) 
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Figure 1 : A pictorial representation of a coarsening process. Blue and pink regions are the growing domains of the two low-temperature equilibrium 
ordered phases. The four panels correspond to snapshots of the system at subsequent times t\ < t2 < £3 < U. 


where g(x,y) [with x = r/L(t w ) and y = L(t)/L(t w )) is a scaling function (and, similarly, s(x) and c(y ) below]. The 
form (0 contains, as special cases, the scaling S (r, t) = s(x) of the equal-time correlation function S (r, t ) = G(r, t , 1) 
and that C(f, t w ) = c(y) of the autocorrelation function C(t, t w ) — G( 0, t, t w ). 

Eq. 0 is, in principle, strictly obeyed only in a quench to a vanishing final temperature Tf = 0. The effect of a 
finite quench temperature is to modify the behavior 0 of correlations at small distances [i.e. r < %{T /), where (-'( Tf) 
is the coherence length of the equilibrium state at Tf (which is small except in the critical region)] and small time- 
differences t - t w < £(Tf) z ]. For simplicity, in the following our discussion about scaling properties will always be 
restricted to sufficiently large space and time scales where such corrections can be discarded. A thorough discussion 
of the whole scaling form taken by G at finite Tf is contained in J2]- 

With the possible exception of liquid systems, whose comprehension however turns out to be more difficult (see 
Article by S. Das, S. Roy and J. Midya in this Volume), real systems are almost never homogeneous, due to the 
presence of quenched lattice defects or because of external disturbances as site-dependent applied fields. Experiments 
show that these sources of inhomogeneities deeply perturb the coarsening process and the scaling properties. Indeed, 
although some kind of scaling symmetry comparable to Eq. 0 is sometimes reported (34| for the correlation func¬ 
tions, the asymptotic growth-law is never observed in the form 0 and, instead, a much slower increase of the ordered 
patterns, usually of a logarithmic type |22l 12 33, 34]), occurs. 

In this Article we will mainly review how and why the growth-law and the dynamical scaling symmetry, expressed 
by Eqs. (fim . are modified in an inhomogeneous coarsening system. 


2. Models 

2.1. Homogeneous systems 

As an optimal playground to study coarsening I consider the description of an homogeneous system in terms of 
the Ising model with Hamiltonian 

'H({S 1 }) = -jY j S 1 S„ (3) 

<‘j> 

where [S ,} is a configuration of N spin variables S, = +1 and i = 1,..., N are sites of a regular (Adimensional lattice. 
By homogeneous it is meant that the Hamiltonian [for instance the one 0 ] does not select any particular spacial 
position and all places are a priori equal. Alternatively one can say that, although a single realization of the system 
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can be inhomogeneous (i.e. due to the presence of domains) local observables do not depend on space after the thermal 
averaging. 

A dynamics is introduced by means of a master-equation with transition rates w({Sj} —> {S'}) regulating the 
hopping between configurations. In the following detailed balance with respect to the equilibrium Boltzmann-Gibbs 
measure e /: ' H will be assumed for the w’s, implying that the equilibrium state is spontaneously approached. I will 
also assume, again for simplicity, that the sequential evolution of single spins occurs and that the order parameter is 
not conserved. 

For such a model, it is well known that, after a quench from an high-temperature equilibrium ensemble to a 
temperature Tf = (kn/jf) 1 below the critical one T c , dynamical scaling holds together with Eqs. (1112b for any Tf <T C 
including Tf - 0. 

In the following I will discuss how the most common sources of inhomogeneities, i.e. varying coupling constants, 
external fields and lattice defects can be included in the above model. I will restrict the discussion to the case in which 
such features are quenched, namely their position, strength etc. is not thermalized and do not vary appreciably during 
the coarsening phenomenon. 

2.2. Varying coupling constants 

If the coupling constants are not uniform the Ising Hamiltonian reads 

ms,]) = - £ (4) 

<ij> 

where — Jq + 0 (/ and the 0’s are usually assumed to be uncorrelated random numbers with OijOu = JJ 2 6,.kdk,(- 
Here and in the following an over-bar denotes an average taken over the realization of the disorder. The probability 
distribution of the 0’s must be such that J t] > 0 (otherwise frustration effects could be present which may radically 
spoil the ferromagnetic character of the system with a drastic modification of the equilibrium and non-equilibrium 
properties). This is known as the random-bond Ising model (RBIM). Notice that, since the distribution of the coupling 
constant is usually assumed to be symmetric around Jq, in this model there is a maximum value Jij - J max < 2 Jo of 
the strength of the bonds. 

2.3. External fields 

Under the action of a site-dependent external magnetic field /;, the Ising Hamiltonian reads 

mSi}) = -jY J S i Sj-Y J h iS i . ( 5 ) 

<ij> 

A useful model is the one where h, is an uncorrelated random variable with expectations hfij = /) 2 b l; , extracted from 
a given distribution (the simplest case is a symmetric bimodal distribution h; = ±h). With these specifications one has 
the so-called random-field Ising model (RFIM). 

2.4. Lattice vacancies 

The effect of quenched vacancies can be taken into account by defining a dilution variable Xij = 0,1 on the lattice 
in such a way that the Hamiltonian reads 


msf) = -j YjXuS.Sj. (6) 

<ij> 

If the dilution variable is factorized as Xij = XiXj the system has site-dilution and the model is usually referred to as 
the site-diluted Ising model (SDIM). Another possibility is to have bond-dilution: In this case Xij is defined on the 
couples of neighboring spins. 

In many cases the ^’s can be considered random numbers, usually uncorrelated in space, such that a fraction D of 
the lattice sites - or bonds - are diluted, namely they are associated to Xi = 0 or^, ; = 0 respectively. This models are 
usually denoted as diluted Ising models. 
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For dilutions larger than the limiting value D c = 1 — P c — where P, is the critical percolation probability - a 
spanning network does not exist. This case is not interesting to us because coarsening implies spins to live on a 
connected graph extending throughout the system. In particular, since D c = 0 in one dimension, these models will not 
be considered in d = 1. 

For the description of other systems, like those where spins are defined on the edges of a fractal network, which 
can be either deterministic or random, then’s are taken in order to reproduce the observed topology of such a network. 
I will briefly discuss this issue later. 

Of course, possible modifications of the models introduced so far - i.e. those where the random variables are 
correlated or where more than one source of disorder is present - are possible, but I will not considered them here. 

3. Downhill versus activated coarsening 

In an homogeneous system described by the model of Sec. 12.11 coarsening occurs at any temperature including 
Tf = 0. This implies that the dynamics, whose main mechanism is a displacement of the interfaces which is regulated 
by their curvature (see article by A.A. Nepomnyashchy in this Volume), proceeds without thermal activation. 

Quenched disorder or any other source of inhomogeneity in the system, instead, usually introduces preferred 
positions where interfaces get pinned in local energy minima. The dynamics can then proceed only by means of 
thermal activation. 

A simple argument shows how this can slow down the kinetics with respect to the homogeneous case where Eq. 
© holds. In view of the dynamical scaling property, the typical pinning energy barrier A E at time t is expected to 
depend on the configuration through Lit) alone, A E = f[L(t), K], where K specifies the set of model parameters, as for 
instance K = (J/TJi/T) in the RFIM or K - ( J/T , d) in the SDIM. In the following, among these various parameters 
we will explicitly write only the one, denoted by e, which denotes the strength of the disorder (as h in the RFIM, or 
2T in the RBIM). 

The evolution is slowed down due to the Arrhenius time t esc (/S.E) oc e^ AE needed to escape pinning barriers. 
Making the simplifying assumption that this effect can be taken into account by a simple rescaling of time ( t —> t/t esc ) 
in Eq. © one has 

Lit, e) = flf 1/2 e~i A£ (7) 

(In the presence of quenched disorder, quantities as L{t) or G(r, t, t w ) are defined throughout as averaged both over the 
thermal history and the disorder realizations). The next point is to establish the form of f[L, e], Although this cannot 
be done in general, most of the systems can be divided into three classes [BJ] according to the way A E depends asymp¬ 
totically on Lit), namely if i) A E approaches a constant value or ii) it diverges algebraically or iii) logarithmically with 
Lit). Let us detail below these cases: 

• i) lim/-,,*, f[L, e ] < c(e), where c(e) is usually an increasing function of disorder (with c — 1 for e = 0), but it is 
constant with respect to L. 

This implies that there is an upper limit to the height of barriers. Using the fact that barriers are approximately 
constant in Eq. © one finds the same growth-law of an homogeneous system but with a smaller pre-factor, of 
order ae~P' if \ which depends strongly on temperature. 

• ii) lim/^oo f[L, e] = b 1 (e)Z/, where (// > 0 is an exponent and hie) is another constant. In this case barriers 
grow algebraically with the domain’s size. This leads to a logarithmic growth 

L(t) ^ [tyT 1 (8) 


for large times. 

• iii) lini/^oo f[L, e] = zie) In L, with zie) a constant. Plugging into Eq. © one finds the large time behavior 

Lit ) oc t ll( , (9) 

with an exponent t) — 2 + fiz which depends on temperature and, possibly, on other system’s parameters. 
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3.1. Role of the lower critical dimension 

The models introduced above [including the usual Ising model ([3} | have a finite critical temperature only in spatial 
dimension larger than the lower critical one d /,. If d < d[ a non-interrupted coarsening process can be observed only 
when quenching to Tf — 0 (in an infinite system), because for any Tf > 0 the ordered phase is not sustained. 

Despite this, by quenching to a very low temperature T /, if the dynamics does not require activation one usually 
does observes coarsening up to a certain time t eq {Tf) such that L(t eq ) — f{T /), where LiTf) is the final equilibrium 
coherence length. For instance, this is what happens in the usual Ising model in d — dr — 1 where f(Tf) - e 2 ^ /J . At 
low temperatures t eq (Tf) grows very large, and for t <sc t eq (T f) the kinetics is analogous to the one at Tf - 0. 

When homogeneity is spoiled the dynamics right at T - 0 is frozen due to the energetic barriers discussed above 
which pin the system in metastable states. The way out is to quench to a temperature high enough to drive activated 
coarsening, but low enough to inhibit the nucleation of equilibrium fluctuations. As for the clean case, as long as L(t ) 
is smaller than the equilibrium correlation length £(T'/), one observes the same coarsening behavior as in the Tf = 0 
quench. 

These considerations apply to all the one-dimensional cases that will be discussed in Section [5] as well to other 
systems, that will be considered in the following, for which T c = 0, as for instance the d = 2 RFIM or some models 
defined on fractal lattices. 

3.2. Crossover structure 

The growth-laws i),ii),iii) of Section [3] are asymptotic behaviors of Lit). Usually these laws are expected when 
L(t ) exceeds a certain characteristic length Aie) associated to the presence of the disorder (or of other sources on 
inhomogeneities). In order to understand the meaning of A(e) by means of a specific example, let us consider the 
SDIM of Sec. 12.41 For D sufficiently small, vacancies are isolated with a typical inter-distance A(D) — I) l/rf . As 
long as L(t) «: A(D ) the domains do not feel the presence of the dilution and the system behaves as an homogeneous 
one (this will be further discussed in Section [7~2T i. Hence, whatever the asymptotic behavior of this model is [either in 
class i), ii) or iii) of the previous Section© this can only show-up for sufficiently large times such that Lit) » A(D). 
One can show that the crossover phenomenon at Lit) - /1(e), where the behavior of the system changes, is a common 
feature, not restricted to the previous example of the SDIM, as we will discuss further. 

On the other hand, the behavior of the systems in the pre-asymptotic regime L(t) <s Aie) is generally model- 
dependent: There are cases where, as in the SDIM discussed above or in models with a fraction of spins frozen in a 
random configuration 0], inhomogeneities are initially ineffective and one recovers the growth-law © of the clean 
case, and others where they produce a different behavior. 

As an example of this different behavior, let us consider the RBIM. In this case there are random couplings on any 
bond, and - if Tf is sufficiently low - this influences the process already at the level of the fast process of single spin 
flips. Due to this, in the pre-asymptotic regime Lit) <s AifJ) the disordered bonds alter the behavior of the system 
with respect to the clean case, thus spoiling Eq. ©, at variance with the SDIM [this will be further discussed in 
Sections l.5. 1 17.31 . However, although determined by the disorder, this early behavior is not the asymptotic one, which 
sets in only when Lit ) is sufficiently large to allow for the simple scaling argument of the previous section[3]to become 
valid. Translating sufficiently large in the condition Lit) » AifJ) provides a definition of Aiff) in this case (although 
physically not as transparent as the one above for the SDIM). 

I conclude this section by mentioning that the slowing down of the growth-law discussed above is not restricted 
to those systems that - when homogeneous - support a non-activated dynamics. The simplest example is provided by 
systems with a conserved order-parameter. In the homogeneous case coarsening proceeds asymptotically by means 
of the activated evaporation-condensation process. However, the presence of disorder [|5t] or other inhomogeneities 
usually introduces further energetic barriers slowing even more the kinetics. 

4. Scaling form 

Once the existence of a crossover length A(e) next to the ordering length Lit) is established, the next step is to 
ascertain if a scaling symmetry is induced by the presence of such lengths and, in case, how the form © can be 
generalized. 
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A straightforward generalization of the form ([2} keeping into account the crossover length is the following 


G(r, t w , t, e) - Q 


Uf) Me) 


L(tw) L(t w ) L(t w ) 


( 10 ) 


and similarly for other observable quantities. This amounts to the usual scaling prescription: To each parameter 
entering G, here r,t w ,t,e, a characteristic length is associated, r, L(t w ), L(t), A(e) respectively, and these quantities 
enters as adimensional ratios in a scaling function Q. 

A debated question regards the dependence of the scaling functions, like @{x, y, z), on the strength of the noise, 
which occurs through the third entry z = A(e)/L(t w ). According to an hypothesis, the so-called superuniversality |@], 
the effect of disorder can be fully accounted for by the slower growth of L(t). This implies that disorder does not 
alter all other properties of the coarsening process, among which - for instance - the geometry of the growing pattern, 
which is encoded in the form of Q in Eq. (fTOl) . This would mean that no extra dependence on e is left in the scaling 
functions and, for the correlations considered above, one should have G(r, t w , /, e) = Q | jjfj, l] = 8 \uj~pj > ^Tj]’ 
where g(x) is the same as that in 0 for the clean system. 

Presently the superuniversality hypotheses has received confirmation in some specific cases 0SSE, HUH, 
but it has also been shown not to be obeyed in other system 113, 14, f5, TfJ 17], as we will discuss below. 


5. Coarsening in one dimension 

The effect of inhomogeneities on the coarsening process can be more easily understood in one-dimensional sys¬ 
tems. In this case the task is 

5.1. RBIM 

Let us suppose to have ites i and i + 1 (namely 

S jS i +1 = -1), as pictorially s □ 



Figure 2: An interface initially located between sites i and i + 1 (plotted in a continuous blue line) moves to a new position r sites away (dashed 
orange line). 

Let us suppose that at a later time the interface moves a distance r away from the original position, namely 
between sites i + r and i + r + 1 (dashed orange line in figure®. The energy change in this process is A E^ +r = 
2[Ji +r j +r+ i - J u+ 1 ]. Since the coupling constants are random numbers this form of A Ejj +r induces local minima and 
maxima in the potential landscape felt by an interface, namely the energy cost of an interface as a function of its 
position, which is sketched in figure® 

Recalling the discussion at the end of Sec. 12.21 there exists a maximum value c(e ) = 4J max for the height of the 
barriers and hence the RBIM is an example of class i) of Sec. [3] Therefore Eq. (Q]l holds for long times (but with 
t < t eq ) at low Tf. 

Clearly, interfaces must travel a certain distance - or, equivalently, domains must grow up to a certain crossover 
size L(t ) = A - in order to sweep a region containing the barrier of maximal height c. According to the argument 
of Sec. 13.21 pinning is effective from the very beginning of the process but the difference with the asymptotic stage 
is that, for Lit) < A , A E may depend on Lit). Accordingly the growth-law of the homogeneous system <jT]) is not 
observed pre-asymptotically and, instead, one finds a ‘J/T -dependent power-law, as in Eq. ® [130. 

Regarding the superuniversality hypotheses it has been convincingly shown not to hold in 111311 . 
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Figure 3: The energy of a single interface is sketched against its displacement r for the RBIM (upper panel) and the RFIM (lower panel) in d = 1 
(we set E = 0 for r = 0). For the RBIM the distance A where the maximum energy bafrrier 4 J max is felt is also indicated. For the RFIM the typical 
energy barrier increases as r 1 / 2 . 


5.2. RFIM 

In this case one has 

i+r 

AE u+r = -Y J h j ~hr l/2 , (11) 

j=i 


where in the last passage we have assumed a bimodal distribution h, = ±h of the random field (but other possible 
choices do not make significative differences) and we have used the central limit theorem. Also in this case there are 
local minima and maxima but the potential landscape where the point defect moves increases with r , as it is pictorially 
shown in figure [3] Indeed, as Eq. (ITU shows, the energy of the interface which has moved to a site at distance r 
can be though itself as the position of a Brownian walker after r steps: the energy landscape is therefore of the Sinai 
type 0. As a consequence, the pinning energy is not bounded as in the RBIM and increases algebraically with the 
typical distance x(t) traveled by the interface, as indicated by Eq. ( fill ). Since it can be shown that also here, as in the 
clean case, the average size of domains L(t) grows proportional to x(f) 0, one concludes that the Id RFIM belongs 
to class class ii) of Sec. [3] One indeed has L(t) oc (In t ) 2 in the large time domain |J4. J_5|. 

Notice that monotonously increasing barriers means that they are bound to prevail over the thermal energy scale 
ksTf if L(t) is large enough, namely for L(t ) » A, where A(h/Tf) ~ (h/T 2 is the crossover length 0 that can be 
obtained balancing the two energetic contributions. 

At early times L(t) « A the Sinai potential is washed out by thermal fluctuations and one recovers the growth-law 
© characteristic of the non-disordered system. The Id RFIM is therefore another example of a crossover with an 
initial homogeneous-like behavior. 

Also in this model superuniversality is not obeyed | (4 15], 


6. Coarsening on fractal substrates 

Coarsening occurring on fractal substrates that can be embedded in an Euclidean space can be still described by a 
model with dilution, like the ones of Section [2~4l I see Eq. ©], with suitable prescriptions for then's. In the following 
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I will focus on the case with site-dilution. The problem has been studied numerically |lj| 20] both for spins defined 
on random fractal networks, like the percolation cluster, and for deterministic ones, such as the Sierpinski fractal and 
others. These deterministic structures are built recursively upon replicating a native geometry - e.g. a triangle for the 
Sierpinski gasket - of a certain size (. 

After the quench, the typical signatures of the coarsening phenomenon are observed, with the growth of an 
ordering-length and dynamical scaling. Correlation functions can be cast in scaling form, as in Eq. <0 : The scaling 
function is different from the homogeneous case and depend on the fractal considered. Regarding the growth-law, 
one observes either a logarithmic increase of Lit), as for a system in class ii) [Sec. [3j, or a temperature-dependent 
power growth-law, as for class iii). This shows that some topological aspect of the substrate - namely the network on 
whose edges the spin variables are defined - is relevant in determining the growth-law and other universal properties. 
However this aspect cannot be traced back to usual indicators such as the fractal or the spectral dimension. 

Despite that a precise identification of this topological feature is currently missing, in |]19(] the conjecture was 
proposed that the same property might determine both the existence of the equilibrium phase-transition and the 
non-equilibrium asymptotic growth-law. Specifically, systems which do not sustain a (finite-temperature) para- 
ferromagnetic transition (examples are the percolation network and the Sierpinski gasket) should belong to class 
iii) and the others (i.e. the Sierpinski carpet) to class ii). 

Concerning the pre-asymptotic behavior, let us stress that, due to the recursive construction of the deterministic 
fractals discussed above, the fractal properties are limited to distances larger than t, while for shorter distances the 
structure is compact. This fixes the crossover length A = l. For Lit) « t one recovers the characteristic growth-law 
0 of a homogeneous substrate. 


7. Coarsening in d > 1 

Understanding coarsening in inhomogeneous systems in dimension larger than one is much harder than in d = 1. 
On the one hand this is due to the fact that analytical approaches are much more difficult, on the other hand because 
physical intuition is less straightforward when interfaces are lines or surfaces that in the presence of the pinning 
centers can bend and stretch. Moreover, numerical simulations - crucial in the absence of analytical tools - are very 
demanding due to the slow growth of L(t). 

Despite many similarities with the one-dimensional models, the pattern of behaviors exhibited in higher dimen¬ 
sions is different and in some cases quite enriched with respect to the Id case, as we will discuss below. 


7.1. RFIM 

There is a general consensus on the fact that the RFIM is characterized by a logarithmic growth both in d = 2 
121, 17, 5 ] and in d = 3 [71, 13), although the value of the exponent i/j defined in Sectioned) cannot be precisely 
determined numerically, due to the very slow increase of Lit). Recalling also the behavior in d = 1 this shows that the 
barriers encountered at late times have a similar nature in all the spatial dimensions considered. Experiments on real 
systems (223 confirm the logarithmic growth. 

The behavior of the model in the pre-asymptotic stage is, instead, more rich. In 13] a pure-like early-stage with Eq. 
(|T} was found. However in |T7] it was shown that, by taking the small T / limit with fixed h/T , a different power-law, 
with an h/Tf dependent exponent C is observed, as in Eq. ([9]». A similar behavior is observed in [23]. 

These apparently contrasting observations can be reconciled in a double crossover scenario. With the two inde¬ 
pendent parameters e = h/T t and t - Tf/J of the model one can build two characteristic lengths Aie) and A(r). Then, 
in principle, one may observe a couple of crossovers as Lit) meets these two lengths, with different behaviors in the 
three time sectors separated by them [furthermore, Eq. Q should be upgraded to contain the extra dependence of the 
scaling function g on A(r)/ L(t w )\. 

In the low-7'/ limit considered in IB where h/Tf is kept fixed, A is finite and r - 0, implying also A - 0. This 
means that there is a single crossover at L(t ) - A from a pre-asymptotic to a late regime where the barriers seeded by 
the disorder act differently as to produce a change from an initial homogeneous-like growth [Eq. 0], to the disorder- 
dependent power-law of Eq. (|9]». In the opposite situation of a large T / (actually as large as possible to keep the 
ordering of the system) A is rather large and A is so small that the crossover associated to it, happening at very early 
times, is not observable. One has again a single crossover, when Lit) - A, but now it occurs similarly to the Id case. 
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namely from a pre-asymptotic behavior where barriers are tamed by thermal fluctuations, when one recovers Eq. (JTJ. 
to the asymptotic logarithmic behavior. 

In the articles cited above, the authors observe alternatively one or the other crossover, and this is possibly due to 
the choice of the model parameters. If the double crossover interpretation is correct it should in principle be possible, 
by tuning the parameters opportunely, to observe the two crossovers in a single quench history. To the best of my 
knowledge, however, this has not yet been reported. 

There is no consensus on the validity of superuniversality in this model. Indeed in fill it was found to hold, 
while it was clearly shown not to be obeyed in & However it must be recalled that the evidences in favor and 
against superuniversality are often obtained with different parameter settings. For instance, the two situations with 
Tf — 0 and Tf finite are considered where, as discussed above, also a different crossover pattern is observed. This 
issue, therefore, is currently not well understood. 


7.2. SDIM 

Numerical simulations of this system in d — 2 have been interpreted both in terms of a power-law increase of L(t) 
19. 24], or as a logarithmic growth-law 125], A reconciliation between these interpretation is provided in [26] in terms 
again of a rich crossover interplay, how it is briefly explained below. 

Coarsening in this system is deeply related to the geometry of the diluted lattice. Let us start, therefore, with a 
brief review of the geometrical properties of the substrate as dilution is varied from high to low values. According 
to percolation theory 127], at D - D c a spanning fractal cluster of non-diluted sites is present, to which an infinite 
percolative coherence length E(Z) C ) is associated (not to be confused with the coherence length of the Ising model 
defined on top of the substrate). This situation is pictorially represented on the far right of figure [4] 



Figure 4: A pictorial representation of the substrate (plum boxes below) as the density D of vacancies is changed (according to the green axis in 
the middle). The behavior of the two lengths A and 3 is shown in orange. 

As D is slightly lowered below D c (i.e. moving a bit to the left in figure [4j, the infinite cluster is fractal over 
distances up to 3( D) ~ (l) c - DY V , where v = 4/3 is a critical exponent, and it becomes compact over larger 
distances. Besides, finite clusters are also present. Upon lowering further D , S keep decreasing until, for a certain 
value D* of D, it becomes comparable to the lattice spacing. The behavior of E( D) is sketched as an orange line in the 
upper part of figure[4][i.e. for D > D*]. For D < IT, the infinite cluster is compact over all length scales and there are 
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no finite droplets. There only remain isolated vacancies and, as already discussed in Sec. 13.21 their typical distance 
defines another length A(D) which keeps increasing upon lowering D up to A(D = 0) = oo, as shown in figure [4] for 
D < D*. 

In the following we show how the conjecture put forth in [191], which was reviewed in Sec. [6] may help in 
understanding the kinetics of this model. Exactly at D c , coarsening occurs on a percolation cluster. Since there 
is no ferromagnetic transition on this network, that is to say that T, - 0, a power-law growth as in Eq. © is 
expected ~~. ' ' ‘ ' “ ~ Q " ‘ ‘ 3int 



Figure 5: A pictorial representation of the behavior of L(t) for different choices of D, on a double logarithmic scale. 


D = D c , the large-scale properties of the network provide an equilibrium phase-transition at a critical temperature 
T C (D ) > 0: In this case, therefore, one expects an asymptotic logarithmic growth-law. This is recognized in figure 
[5]by the asymptotic downward bending of the curves with 0 < I) < I), for large times. Clearly, right at I) = 0 the 
homogeneous case is recovered with the usual law (Q}, which is shown in red in figure 0 

This pattern of behaviors is indeed observed in the numerical study 112611 of the model in the small T /-limit. This 
supports the conjecture discussed in Section[6] 

Not only the asymptotic laws, but also the early stage behaviors can be understood. Let us start from the large 
dilution regime D < D c . In this case the properties of the spanning cluster are the same as those at D, up the length 
S(D). Hence a crossover phenomenon is expected from an early stage regulated by Eq. <0 for L(t) <s E(D) to the 
asymptotic logarithmic growth when L(t) » E( D), as it is pictorially represented with a dotted-orange line in figure 
0 This shows that the growth-law may be interpreted as algebraic or logarithmic depending on the choice of D and 
on the time-scale of the observations and might explain why different findings are reported in the literature. 

Let us turn now to consider the small dilution regime, which behaves differently. As discussed in Section l3T2l a 
crossover occurs here between an early homogeneous-like behavior © and a late logarithmic growth when L(t) meets 
/1(D) ~ D -1 ^. This is shown in figure0with a dotted-violet line. 

This whole pattern of crossovers is actually observed in [26] suggesting - once again - that an interpretation in 
terms of the scaling symmetry and the related crossover structure is a groundbreaking tool. 

Notice also that, although two lengths are present also in the SDIM, the crossover mechanism is very different 
from the one - denoted as double crossover in Section [7TTI - invoked for the RFIM. Indeed, in that case the two lengths 
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A(e) and A(r) are built from two independent parameters of the model. This allows one, in principle, to have a double 
crossover, namely three consecutive regimes. 

For the SDIM in the small-7 'j limit considered in [26], both 3 and A depend on the same parameter D and this 
occur in such a way that either A - 0 (for D > D*) or A =* 0 (for D < D*). Hence, for any choice of 0 < D < D c , only 
one between the two lengths A and A can play a role, and a single crossover between an early and a late stage can be 
observed, as it is clarified in figure [3 

Finally, let us consider the issue of superuniversality. In Ref. J25| [24] the authors succeed in collapsing the 
curves for correlation functions using the simple scaling (0. In [25], this is interpreted as a possible confirmation 
of the superuniversality hypotheses. However in [26] a study with a more ample choice of values of D and different 
timescales clearly shows that (0 is not sufficient to collapse all the curve and, instead, a scaling form with an extra 
argument as Eq. lITOl) is needed, suggesting that superuniversality is violated also in this model. 


7.3. RBIM 


For the RBIM, Huse and Henley (HH) [28]], elaborating on previous ideas due to J. Villain [29], predicted the 
logarithmic law 0 with ijj — 1/4. Their argument, in d - 2, can be summarized as follows: According to roughening 
theory a piece of interface of linear size r in equilibrium deviates from flatness of a typical quantity w(r) oc A', where 
the roughening exponent ( r is a universal quantity. If disorder is present, the pinning energy of the interfaces scales as 
E p ~ r x , where \r is another exponent. 

On the other hand, simply due to geometry, moving a distance r along the surface of a coarsening curved domain 
one also deviates from flatness of a certain quantity a(r,R), where R is the radius of curvature. Hence, on scales r up 
to r c - such that a(r, R) - w(r) - the domain’s wall is basically equilibrated. A simple calculation shows that air, R) is 
at least of order a = r 2 /R (this can be very easily proved for a circular shape in the small r limit) and hence one has 

r c * Ri/C-M' 

The next smallest scale to evolve towards equilibrium is precisely r c and, recalling the discussion above, the 
pinning energy associated to this is E p ~ rf. r - R' ! ', with if/ — XrtO- ~ tr)- For random bond disorder in d = 2 it is 
known [28, 30] that = 2/3 and Xr = 1/3, leading to <// = 1/4. Assuming scaling by the identification R ~ Lit) and 


plugging this result into ii) of Section [3]one arrives at Eq. ® with i/y = 1/4. 

The correctness of the HH prediction has remained controversial for quite a long time since numerical results 


for the d — 2 RBIM showed a steady algebraic growth with a disorder-dependent dynamical exponent [81 [31]] as in 


Eq. 0. This would imply an asymptotic logarithmic increase of the pinning barriers at variance with the prediction 
of HH. Recently, more extensive numerical simulations d have produced sufficient numerical evidence for the 
existence, after a long-lasting algebraic regime, of a crossover to an asymptotic logarithmic growth. This supports the 
HH scenario, although a precise determination of the ijj exponent is, also in this case, not possible. Experiments on 
real two-dimensional systems also find 132, 33l 34] a logarithmic growth and the HH exponent i/j = 4 was reported in 

[HH. 

The very nature of the long-lasting pre-asymptotic regime with power-law increase of the domain’s size, which 
is observed also in some experiments 034ll . is not presently clarified. In view of the discussion of Section Q] it could 
be due to a genuine regime where the scaling of barriers is only logarithmic - according to iii) - different therefore 
from the asymptotic behavior described by HH. This is not necessarily the case, since even within the HH scenario 
long-lasting pre-asymptotic corrections might yield an effective algebraic growth-law at intermediate times, as noticed 
in different contexts in [35]]. 

Interestingly enough, the algebraic growth-law can be made truly asymptotic by considering a symmetric bimodal 
distribution J,j = Jo ± Jo of the coupling constants, as it was shown numerically in [ 36]. This choice amounts to have a 
fraction 1 /2 of bonds , = 0 and an equal quantity of J t j = 2 Jo. This realizes a percolation cluster of bonds, since the 
bond percolation threshold is D c = 1/2. According to the conjecture mentioned in Sec. [3 one expects an algebraic 
increase of L(t) as in Eq. 0, as indeed it is observed. 

Regarding super universality, also in this model the situation is controversial: while it was observed to hold true in 
Sill, clear violations have been observed in ifTfih . Some experimental evidences in favor of a superuniversal scaling 
have been reported in [34]. 
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8. Conclusions 


In this Article I briefly reviewed the coarsening phenomena occurring in simple model-systems without space- 
translation invariance. Following in the footsteps of early studies of phase-ordering in homogeneous ferromagnets, 
most efforts in this field have been devoted to the determination of the growth-law and of the dynamical scaling 
properties. Apart from some instances of one-dimensional systems, analytical approaches are scarce and the current 
understanding of the problem rests mostly upon heuristic arguments and slowly converging numerical simulations. 
This fact, together with the observation of a broad spectrum of distinct behaviors, make the problem far from being 
settled. 

Our poor understanding hinders the development of a general framework for growth-kinetics in disordered and 
inhomogeneous ferromagnetic systems, which is a subject of noteworthy thechnological relevance. This, in turn, en¬ 
cumbers the discussion on the very nature of more complex disordered systems, such as spin-glasses, which according 
to a debated interpretation 0 could also be described as a disordered ferromagnet. 

As discussed in this paper, despite this state of affairs, the rich (and sometimes apparently contrasting) behavior 
of different systems - or even of the same system in different temporal or parameter regions - can be interpreted in 
terms of an upgraded scaling framework where, next to the typical size L of ordered regions, there is another (or, in 
some cases, others) relevant characteristic length A associated to the inhomogeneities. This foster the overview of 
the phenomenon and could help to set the basis for a general renormalization group scheme where diverse univer¬ 
sal behaviors - e.g. power-law or logarithmic growth-laws - are connected by a crossover pattern. However, the 
development of such a controlled theory is matter for future research. 
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